Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 3 February 2008 (MN WF$t style file v2.2) 



Constraints on the nature of dust particles by infrared 
observations 



Cs. Kiss 1 , P. Abraham 1 , R.J. Laureijs 2 , A. Moor 1 , S.M. Birkmann 3 

i 1 Konkoly Observatory of the Hungarian Academy of Sciences, P.O. Box 67, H-1525 Budapest, Hungary 
i 2 European Space Agency, Astrophysics Division, Keplerlaan 1, 2201AZ Noordwijk, The Netherlands 
3 Max-Planck-Institut fur Astronomie, Konigstuhl 17, D-69117, Heidelberg, Germany 



Accepted: 2006 September 24; Recieved: 2006 September 23; in original form: 2006 April > 



ABSTRACT 

The far-infrared (FIR) emissivity of dust is an important parameter characterizing 
the physical properties of the grains. With the availability of stellar databases and 
far-infrared data from Infrared Space Observatory (ISO) it is possible to compare the 
optical and infrared properties of dust, and derive the far-infrared emissivity with 
respect to the optical extinction. In this paper we present the results of a system- 
atic analysis of the FIR emissivity of interstellar clouds observed with ISOPHOT (the 
photometer onboard ISO) at least at two infrared wavelengths, one close to ^100/im 
and one at 200/^m. We constructed FIR emission maps, determined dust tempera- 
tures, created extinction maps using 2MASS survey data, and calculated far-infrared 
emissivity for each of these clouds. We present the largest homogeneously reduced 
database constructed so far for this purpose. During the data analysis special care 
was taken on possible systematic errors. We find that far-infrared emissivity has a 
clear dependence on temperature. The emissivity is enhanced by a factor of usually 
less than 2 in the low dust temperature regime of 12K<Td<14K. This result suggests 
larger grain sizes in those regions. However, the emissivity increase of typically below 
2 restricts the possible grain growth processes to ice-mantle formation and coagulation 
of silicate grains, and excludes the coagulation of carbonaceous particles on the scales 
of the regions we investigated. In the temperature range 14 K < T<j < 16 K a systematic 
decrease of emissivity is observed with respect to the values of the diffuse interstellar 
matter. Possible scenarios for this behaviour are discussed in the paper. 
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1 INTRODUCTION 

The far-infrared (FIR) emissivity of interstellar grains pre- 
dominantly emitting at wavelengths in excess of 100 fim can 
now be determined in many molecular and moderate density 
regions thanks to the increasing availability of far-infrared, 
sub-mm, and massive stellar data sets. Several studies have 
reported on the detection of an enhancement of the FIR 
emissivity in the 100-200 /im wavelength range in regions of 
higher column density compared to the dust emissivity in 
the diffuse interstellar medium associated with HI (Bernard 
et al., 1999; Cambresy et al., 2001; Juvela et al., 2002; Step- 
nik et al., 2003; Lehtinen et al., 2004, Rawlings et al., 2005). 
The largest sample of different regions is presented in del 
Burgo et al. (2003) presenting a study of eight translucent 
clouds. 

The observations suggest a trend where the emissivity 
increases with decreasing temperature. The variation is at- 



tributed to a change in grain properties which is expected to 
take place in denser environments. In particular, the increase 
in emissivity is interpreted as a signature of an enhancement 
in grain size. 

However, the trend shows a large scatter which might 
be due to the observational limitations which increase the 
uncertainties. On one hand, star counts statistically probe 
only a limited extinction range and, on the other hand, the 
determination of the grain temperatures in the infrared can 
have large uncertainties. 

In this study we present a large sample based on In- 
frared Space Observatory (ISO) data of cloud regions with 
reliable values of the FIR dust emissivity and optical extinc- 
tion data. We examine in detail the possible observational 
and data processing errors to ensure the reliability of our 
results. A general good agreement was found between our 
results and previous studies of individual regions. Due to 
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the large number of data points we are able to put stronger 
constraints on the emissivity changes with temperature. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Far-infrared maps 

We searched the ISO Archive (Salama 2004) for ISOPHOT 
observations (Lemke et al. 1996) of interstellar clouds match- 
ing the following criteria: (1) the field has been covered at 
least at two far-infrared wavelengths: one at 200 /im (long) 
and another either at 90, 100 or 120 /im (short) in order to 
provide a sufficient wavelength interval for a reliable colour 
temperature calculation; (2) the cloud has to be galactic and 
must have sufficient dynamic range in brightness for corre- 
lation analyses; (3) the map is larger than 5' at least in 
one dimension; and (4) there is no high mass star formation 
going on in the vicinity of the cloud which could signifi- 
cantly change the local interstellar radiation field. In total 
we selected 22 maps which is the largest sample evaluated 
so far for studying far-infrared dust emissivity. All selected 
maps were obtained with the P22 astronomical observing 
template mode (Laureijs et al. 2003). Measurement wave- 
lengths, ISO-id numbers and central positions are listed in 
Table 1. 

The ISOPHOT observations were performed with the 
C100 (43" 5x43'.' 5 sized pixels) and C200 (8975x89'.' 5 sized 
pixels) cameras. The ISOPHOT data were processed with 
the Phot Interactive Analysis software version 10.0 (PIA, 
Gabriel et al., 1997), using standard batch processing and a 
first quartile flatfielding. We followed in detail the processing 
scheme described in del Burgo et al. (2003). The data were 
colour corrected taking into account the dust temperature 
derived from the brightness ratio at the two wavelengths 
(see Sect. 3.1). 

The officially quoted absolute photometric uncertainty 
of the surface brightness in the ISO Legacy Archive is <20- 
25% (Klaas et al., 2003). In the present study we esti- 
mated independently this uncertainty in two ways. We com- 
pared the ISOPHOT and COBE/DIRBE background sur- 
face brightness values, interpolated to the ISOPHOT wave- 
lengths (1) for a large sample of mini-map observations 
and (2) for our target fields. This analysis is presented in 
Appendix A. The results of this investigation show, that 
the typical relative deviations of the individual ISOPHOT 
measurements with respect to the COBE/DIRBE values is 
~15%, and that there is no noticable systematic discrepancy 
between the two photometric systems. Throughout this pa- 
per we present our results in the ISOPHOT/PIA 10.0 surface 
brightness photometric system. The effect of a potential im- 
perfect surface brightness calibration is discussed in detail 
in Sect. 4.5, based on results presented in Appendix A. 



2.2 Extinction Data 

The traditional method to derive extinction in interstellar 
clouds is based on variations in stellar density in the sky due 
to the obscuration of dust (Wolf 1923). Star counts obtained 
by placing a regular grid on the target field are converted to 
B or V-band extinction using statistical methods; the zero 



extinction level is obtained by comparison with a nearby 
extinction free reference field. 

Cambresy et al. (1997) replaced the classical regular 
grid with an adaptive one; in this method the gridsize is 
adjusted to include a fixed number of stars, therefore it can 
still provide extinction estimates for high density regions. 
However, in these cases the derived Av values are averages 
over the enlarged area of the adaptive cell. 

Due to the requirement of sufficient count of stars for 
reliable statistics, the minimum resolution is limited to ar- 
cminute scales, and low dust column densities. The max- 
imum extinction in the visual is at best 5 mag and cannot 
be improved significantly by deep dedicated observations. In 
that respect, online catalogues with optical stellar data are 
well suited for the method. 

In recent years the application of near-infrared redden- 
ing of individual stars has become widely used due to the 
general availability of J, H and K-band measurements. Ex- 
tinction mapping methods using NIR reddening often com- 
bine the individual reddening values in some statistical man- 
ner in order to minimize the effect of individual line of sights. 
The NICE and NICER colour excess methods (Lada et al. 
1994, Lombardi & Alves 2001) proved to produce good qual- 
ity maps in many applications. Although NIR star counts 
are more reliable for high extinction values, in the low to in- 
termediate extinction range (Av < 15 mag) the NIR colour 
excess methods are superior over NIR star counts, since 
they are less affected by the presence of foreground stars 
(Cambresy et al. 2002). 

The 2MASS Point Source Catalogue (Cutri et al. 2003) 
is a powerful source of NIR data and can be used to obtain 
NIR colour excess / extinction as was done by many authors 
in recent years. 

A detailed comparison of visual extinctions obtained 
by optical star count and NIR colour excess is presented in 
Kiss et al. (2006) 1 . In this paper extinction maps are derived 
using star counts of USNO data (Monet et al., 1998, 2003) 
and using near-infrared (J, H and Ks-band) reddening of 
2MASS data. Their fig. 2 shows a good correlation of the two 
datasets and the slope of the scatter plot is very near to unity 
for low extinctions. A flattening of the USNO extinction 
occurs at 3. ..3. 5 mag, due to the already insufficient count 
statistics at these extinction levels. 

The near-infrared colour excess methods are less af- 
fected by uncertainties due to foreground stars and give re- 
liable extinction values in denser regions up to ~15mag. 
Therefore we decided to apply the NICER method on the 
2MASS Point Source Catalogue to obtain extinction maps in 
this paper. Following the findings by Cambresy et al. (2002) 
we set Ay" 1 = 8 mag for the liming magnitude of 2MASS 
data, assuming a 10% of foreground stars in the target field. 
Above this limit we did not consider the extinction data to 
be reliable. These data points were excluded from the fur- 
ther analysis. Due to the low number of data points above 
this limit and outlier resistant routines we used in the scat- 
ter plot analysis, difference between the effect of a cut in 
Av or of a cut perpendicular to the slope at Ay m = 8 mag is 
negligible. 



1 Note, that the first author of the cited paper is not the first 
author of the present paper 
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Figure 1. Example (a) of a short versus long wavelength surface 
brightness scatter plot and (b) a visual extinction versus 200 (im 
scatter plot, both in the G100.0+14.8 field. On both panels the 
dashed lines show the fitted linears. 



3 RESULTS 

3.1 Dust temperature 

We derived colour temperatures Td from the slope in the 
l\ 1 versus Ia 2 scatter plot of a specific sky region, where I\ ± 
and I\ 2 are the surface brightness values at the wavelengths 
Ai and A2, respectively. An example in shown in Fig. la. 
Maps taken with the C100 camera (90 and 100 /im) were 
smoothed to the resolution of the 200 (im observation. The 
method of using scatter plots to determine Td is described 
in del Burgo et al. (2003). It has the main advantage that 
the slope is insensitive to surface brightness offsets due to 



e.g. zodiacal light, extragalactic background, calibrational 
zero point, etc. The conversion of the AIaj/A^oo slopes 
to T d was done via high resolution AI^/AI^oo- T d tables, 
assuming a i/B„(T) spectral energy distribution. The tables 
also accounted for the correct colour correction. We adopted 
/? = 2. 

In Sect. 4.3 we check the effect of a different emissivity 
law. In the following we adopt Td as an approximation of 
the physical temperature of the dust grains. 

Del Burgo et al. (2003) found that the scatter plots 
of some fields can consist of more than one linear section, 
indicating multiple dust temperatures within the field. In- 
spection of the maps of these fields indicated that the differ- 
ent temperature components came from interlocked regions 
which could be separated by Av isocontours. The resulting 
Iai versus I200 scatter plots of the subfields were sufficiently 
linear. 

The uncertainty in the temperature Td can be at- 
tributed to, firstly, the error in the slope 5(AIaj /AI200), and, 
secondly, the calibration uncertainty in the surface bright- 
ness values. The first component was computed from the 
formal uncertainty of the linear fits. The second component 
was evaluated by using surface brightness uncertainty val- 
ues given in Sect. 2.1. The resulting temperature and uncer- 
tainty values are typically in the order of 0.7 K depending 
on the colour temperature (Table 1.) 

3.2 Emissivity parameters 

In Table 1 we list the derived values of I200/AV obtained 
from the slopes of the I200 versus Av relationships in those 
regions where Td is constant as delineated by a constant 
1\ 1 vs. I200. The 200/im surface brightness was smoothed 
to the resolution of Av- The ratio I200/AV consists of two 
independent observables and therefore also serves as a good 
diagnostics of the main uncertainties. An example Av vs. 
I200 scatter plot is presented in Fig. lb. 

Using the measured value of Td from the ratio 1200, 
we derive the ratio between the infrared and optical opacity 
T2oo/Av from T2oo/A v =(l20o/A v ) x B 2 oo(T d )~\ where 
B20o(Td) is the Planck function at 200 /im at the dust tem- 
perature T d . 

Altogether, we derived 26 sets of T d , I200/AV and 
T2oo/Av values in the 22 maps, due to multiple tempera- 
ture components in 4 fields. In addition, six fields of the del 
Burgo et al. (2003) sample have been reprocessed following 
our scheme (see Table 1). 

The errors of T200 /Av come from two sources: from the 
error in the determination of I200 /Av in the I200 versus Ay 
scatter plot and from the temperature uncertainty in B„(T). 
Assuming that these two sources are independent, the final 
T2oo/Av errors can be expressed via the partial derivatives 
by I200/AV and B„(T), following standard error propaga- 
tion. Since we use the slope of the I200 versus Av relation 
only, systematic offsets in I200 and Av - like in the temper- 
ature computation - do not play a role. 

To determine the propagation of the error in T d in the 
Planck function, we use the two values T-<5T and T+<5T to 
estimate the upper and the lower error bars in T2oo/Av- In 
most cases the upper and lower error bars are nearly similar, 
therefore we give one average value for the uncertainty in 
Table 1. In the subsequent figures we present these upper 
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and lower error bars individually. The relative errors caused 
by the temperature uncertainties are typically much larger 
than the relative errors from the determination of 1200/Av 
in the scatter plots. 



3.3 Emissivity versus temperature relationships 

In Fig. |3 we present the resulting ratios of 1200/Av and 
T20o/Av as a function of Td for our sample of targets, includ- 
ing the reprocessed del Burgo et al. (2003) sample. In ad- 
dition, we included in Fig. 03 the values obtained by differ- 
ent authors (for references see figure caption). The I200/AV 
ratios of the different studies were made consistent by con- 
verting to A= 200 ;Um assuming a modified black-body with 
a (3 — 2 emissivity law. The T200 / Av values from other au- 
thors, which were originally obtained at longer wavelengths, 
were transformed to 200 zjm as well. A comparison of our 
results with these previous works is presented in Sect. 5.1. 

In Fig. |21 the dashed line marks the value of the dust 
emissivity representative of the DISM. Note, that the DISM 
itself has a unique temperature of ~17.5K. The observed 
values of T2oo/Av exhibit significant deviations from this 
reference value as a function of temperature. Fields above a 
dust temperature of 14 K show lower emissivity value than 
that of the DISM with no obvious trend in their distribution. 
In the colder fields emissivity values are above that of the 
DISM, with a weak indication of a decreasing trend towards 
lower temperatures. 

A few points are above the typical emissivity values 
found for the majority of the fields at low temperature. The 
higher uncertainties in these fields - as derived from the un- 
certainites of the temperature and I200 / Av ratio determina- 
tion - cannot explain the high T2oo/Av values observed. Fur- 
ther potential sources of systematic errors that could lead to 
incorrect emissivity values are discussed in detail in Sect. 4. 



4 UNCERTAINTIES IN THE COMPUTATION 
OF EMISSIVITY PARAMETERS 

The reliability of the final emissivity values derived for our 
target fields may be effected by several sources of errors. 
This is specially true, if a certain source of error does not 
contribute only to the random scatter of the derived values, 
but causes a certain pattern, which may be mistaken for a 
physical change in dust properties. In this section we inves- 
tigate several potential error sources, and try to estimate 
their impact on the final emissivity determination. 

4.1 Derivation of the correct extinction value 

The extinction value Av depends on the wavelength regime 
and the method used to derive it. In Sect. 2.2 we compared 
the results from the USNO star count and the near-infrared 
NICER method and concluded that the latter one is more 
suitable for our study. 

The presence of foreground stars may also have an ef- 
fect on the derived visual extinction values. While in a low 
extinction part of the cloud most stars in the beam are 
background objects, towards the densest peaks the stellar 
sample in the beam may be dominated by foreground stars. 
As a result, the visual extinction in the densest parts can 
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Figure 2. (a) Relationship between I200/AV and the 
dust temperature. Dashed lines represent iso-emissivity 
(e =[t"20o/Av]/[t20o/Av] DISM = constant) curves. The cor- 
responding emissivity values are indicated, (b) Relationship 
between the 200 fim emissivity (T200 / Av ) and the dust colour 
temperature. Black dots with error bars mark the fields inves- 
tigated in this study. Open circles with error bars mark the 
reprocessed del Burgo et al. (2003) target fields. We also present 
recent results published in other papers, marked by gray symbols: 
- filled triangles: Stepnik et al. (2003); — filled, upside-down 
triangle: Pagani et al. (2003); - filled squares: Rawlings et al. 
(2005); — cross: Bianchi et al. (2003); — open triangle: Lehtincn 
et al. (1998); - open diamonds: Lehtinen et al. (2004); — open 
square: Juvela et al. (2002). 



be underestimated, and the emissivity value overestimated. 
Since the foreground objects show low E(B-V) values, we 
minimized the bias they cause by using an outlier resistant 
method to calulate the average reddening value in each cell 
in the NICER method. However, all of our target clouds are 
nearby (closer than a few hundred parsecs) and are located 
at high Galactic latitudes (b>15°), therefore the contribu- 
tion of foreground stars is negligible. 

The extinction value may also depend on the size of the 
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Table 1. Summary table of the basic properties and derived quantities of the fields analysed in this study. The columns of the table 
are: (1) the name of the field (denoting central galactic coordinates); (2) ISO-id (TDT-number, Laureijs et al., 2003) of the short and 
long wavelength measurements. (3—4) right ascension and declination of the field centres (J2000); (5) short and long ISOPHOT filter 
wavelengths used; (6) dust temperature derived from ISOPHOT short- versus long wavclenght (200 (im) surface brightness scatter plots, 
the uncertainties presented here reflect the formal errors of the fits only; (7) fitted slope of the 200 fim surface brightness (I200) versus 
visual extinction (Ay) scatter plot; (8) dust emissivity derived from columns #6 and #7. 



grid and whether it is of fixed or adaptive type grid. Cal- 
culation of extinction/reddening using an adaptive gridsize 
may significantly underestimate the extinction, especially in 
a cloud centre, due to a inhomogeneous distribution of stars 
in the cells. We therefore preferred to use a fixed gridsize in 
our analysis. 

4.2 Dependency of Ry 

The ratio of total over selective extinction Rv is found to be 
significantly higher in dense interstellar clouds than Ry in 
the diffuse ISM (Cardelli et al., 1989; Whittet et al., 2001). 
This change is a strong indicator that the optical proper- 
ties of dust grains differ in dense regions. When we derived 
Ay from Aj by following the method of Lombardi & Alves 
(2001), we adopted Ry = 3.1, which is a value typical for the 
diffuse ISM and may not be representative for the regions 
we investigate. 

Applying the empirical relationship between Ry and 
Aj/Ay as determined by Cardelli et al (1989), we find 
Aj/Ay = 0.282 and 0.334, for Ry = 3.1 and 5.5, respectively, 



implying that Av will decrease by 19% when increasing Ry 
from values representative of the diffuse ISM to dense lines 
of sight. The corresponding I200/AV and T2oo/Ay values will 
increase accordingly by 19%. We infer that the uncertainty 
in the assumption of Rv = 3.1 in the determination of Ay 
has a relatively minor implication on the values of the ob- 
servables presented in this paper. The assumption may in- 
troduce an underestimation of at most 20% in both I200 /Ay 
and r/Ay. 

4.3 Application of a different emissivity law 

Assuming that the exponent of the emissivity law /3 is 
close to 2 in the interstellar medium (Gezari et al. 1973, 
Draine & Lee, 1984, Boulanger et al., 1996), we adopted 
/3=2 to derive T<j from the slope of the short (90, 100 
or 120 fim) versus long (200 /an) wavelength scatter plots 
(Sect. 3.1). Based on balloon-borne submillimetre observa- 
tions of a large number of ISM regions, Dupac et al. (2003) 
found that /3 is not constant but depends on the dust tem- 
perature as /3(T) = (0.4 + 0.008 ■ (T/l K)) _1 . We used this 
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Figure 3. Displacement of t/Av vs. Td data points due to the 
effect of a temperature dependent /3 (Dupac et al., 2003). The 
original emissivity values with (3 = 2 are marked by filled circles. 
The stripes originating from these points indicate the change due 
to a different (3. 
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Figure 4. Displacement of the Td-T2oo/Av data points with the 
correction for the warm dust component of the short wavelength 
surface brightness values. The original data points are marked by 
filled circles and the stripes originating from them indicate the 
change due to the correction for the presence of warm dust. 



relationship to investigate the effect of possible changes in 
(3 on our results. We calculated the new colour tempera- 
ture from the I s hort versus I200 scatter diagrams using /3(T). 
From this temperature we obtained the modified values for 

T20o/Av- 

In Fig.|3]we show the effect of /3{T) on the temperature 
versus emissivity relationship. The /3 dependency tends to 
increase the temperature, but the change is small, only for 
temperatures higher than 15.5 K the variations are notice- 
able, more than 0.3 K but less than 0.5 K. Although the new 
temperature-dependence of (5 cause a systematic change of 
the distribution of the measurements in Fig. |3J the effect is 
small with respect to the measurement uncertainties. 

4.4 Presence of warm dust 

IRAS related studies of the dense and diffuse ISM have 
shown that extended emission at 60 /im can be associated 
with dust in low density regions. In these regions the 100 
/im emission is closely correlated with the 60 fim emission, 
suggesting that they trace the same dust component in the 
ISM. At higher densities the 60 (im emission becomes weak 
indicating the absence of grains giving rise to the 60 ^m 
emission. The remaining 100 /im emission is a signature of 
a "cold" dust component at T <15 K. This observation is 
supported by the close correlation of the cold dust compo- 
nent with 13 CO emission. The surface brightness measured 
at a wavelength of A ~ 100/im may contain emission from 
both the cold and the "warm" (or low density) dust compo- 
nent. The warm component will affect the determination of 
the I v / 1200 colour temperature, such that the temperatures 
will be overestimated due to the short wavelength emission 
component. 

To determine the fraction of the cold component in 
our short wavelength (~100/im) maps we analysed 60 and 



100 fim maps of our regions obtained with IRAS /ISS A be- 
cause in most fields there are no ISOPHOT observations at 
60 fim. The amount of emission from the cold component 
at 100 fim towards a specific sky position can be calculated 
from the two IRAS bands: 

Icoid = Ixoo — 8 ■ leo (1) 

where 9 is the slope of a linear relationship fitted to the 
60 vs. 100/j,m surface brightness correlation diagram in the 
outer parts of the cloud. The relative contribution of the cold 
dust component at 100 fim is estimated from the correlation 
between l co id and the original I100 surface brightness. The 
scatter diagram was fitted by a linear function, resulting in a 
slope X100 , which gives the ratio of the cold dust emission to 
total surface brightness at 100 fim. We find typical values of 
0.7 < X100 < 0.8. We used X100 to determine the cold dust 
emission contribution in the ISOPHOT short wavelength 
bands. Since the ISOPHOT measurements were not always 
taken at 100 (jm, we converted X100 to X\ assuming that the 
spectral energy distribution of the total surface brightness 
can be approximated by two modified black-bodies (/3 = 2) 
with two different temperatures, T marm and Tcoid- We fixed 
Twarm to 17.5K (Lagache et al., 1998) and set T co id to the 
temperature obtained from the original surface brightness 
scatter plots. Our analysis shows that X\/Xioo is close to 1 
and the ratio is not very sensitive to T co id- These X\ values 
were applied to the short wavelength surface brightness val- 
ues at 90, 100 or 120 fim, to compute the corrected T2oo/Av 
values. The results are presented in Fig.0] Fig.0]shows that 
the presence of warm dust tends to overestimate the tem- 
perature by 0.5-1 K which is significant. As a result, due to 
the lower temperature of the cold dust and since /200/Av is 
hardly affected, correction for warm dust emission increases 
warm dust increases T2oo/Av by at most 30%. 
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Figure 5. Effect of calibration errors in the surface brightness cal- 
ibration, (a) Statistical errors in the short and long wavelength 
surface brightness (see the text for details). The four lines orig- 
inating from each black dot (the original T2oo/Av data points) 
span a banana-shape area in the figure, corresponding to 9% cal- 
ibration errors in both the short and long wavelength ISOPHOT 
surface brightness values, as found in Appendix A. (b) Displace- 
ment of the T200 /Ay values (stripes) relative to the original ones 
(black dots, as listed in Table 1), when the correlation coefficient 
between the ISOPHOT and DIRBE surface brightness calibration 
(see Appendix A) are applied. 



4.5 Effect of surface brightness calibration errors 

Errors in the ISOPHOT surface brightness calibration may 
affect the derived dust temperatures and the relationships 
presented in Fig.|5] Here we first discuss the general effect of 
surface brightness calibration errors applying a simple model 
in Sect. 4.5.1. 

Actual calibrational errors in this photometric system 
can only be unraveled by a comparison with a standard 
system. Due to its measurement design the DIRBE instru- 
ment on-board the COBE satellite was able to perform accu- 



rate absolute photometric observations (Hauser et al. 1998a) 
and serves now as the standard system for infrared sky 
brightness observations. A detailed comparison of the sur- 
face brightness photometric systems of ISOPHOT and that 
of COBE/DIRBE is presented in Appendix A. In Sect. 4.5.2 
we apply the ISOPHOT-DIRBE surface brightness trans- 
formation equations found in Appendix A, and discuss their 
effect on our results. 



4-5.1 General effect of calibration errors 

In order to test the effect of the surface brightness cali- 
bration errors on the derived emissivities, in a simple ap- 
proach we assumed a linear relationship between the true 
sky brightness and the surface brightness measured by a 
specific instrument/filter, i.e. the transformation equation 
between the two can be written as: 



= G A xJf y + C 



(2) 



where I™ oas is the measured sky brightness, I^ rue is the 
real sky brightness, and G (gain) and C (offset) are the co- 
efficients describing the transformation. 

Our analysis method, which is based on slope fitting in 
scatter plots (see Sect. 2), is not sensitive to errors due to 
constant offsets in the surface brightness. Multiplicative er- 
rors, however, change both the colour temperature and the 
emissivity parameters. Such errors may come from the pres- 
ence of a gain uncertainty due to incorrect surface brightness 
calibration. We investigated the impact of this effect on our 
emissivity values, using a simple model. We assumed that 
I s m ho a r s t =G short -I^ rt and l?oo as = G 2 ood s 2 k ou , where I^ rt and 
Ijqq are the true surface brightness values of the sky at our 
short wavelength (either 90, 100 or 120 fim) and at 200 fim, 
respectively, while Ishort and l2uo as are the values affected by 
calibration errors. G s hort and G200 are the detector gains for 
the short wavelength and for the 200 fim filters, respectively. 
For both Gshort and G200 we assumed two extrema, 0.91 and 
1.09, correspondig to an uncertainty value of 9%, as found 
for the ISOPHOT filters in Appendix A. 

The effect of the four gain combinations for the de- 
rived emissivity values and dust temperatures is presented 
in Fig. 5a. The shorter, fully vertical lines originated from 
the black dots correspond to the cases, where both Gshort 
and G200 have the same value (either the same 0.91 or 1.09). 
In these cases the changes in the T2oo/Av values are 9%. 

The longer, tilted lines correspond to the cases, when 
Gshort and G200 have opposite values. Then the amount of 
change in T2oo/Av and Td depends noticeably on the origi- 
nal temperature, but in general is 25. ..30%. This is a major 
effect, but the T2oo/Av values remain in the domain, that 
their relative deviations from the emissivity of the diffuse 
interstellar matter is still significant. 

The difference between the original and recalcultated 
dust temperatures are in the order of ~0.7K. As discussed 
in Sect. 3.1 these calibration errors are the main sources of 
uncertainty in the temperature determination. 

The scenario presented here is indeed a 'worst-case', i.e. 
the relative effect of the calibration errors at the two wave- 
lengths are the strongest. The likely effect is smaller and 
would result in a final T2oo/Av uncertainty of 10. ..20%, de- 
pending on the original colour temperature. Here we also 
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did not take into account the systematic calibrational dif- 
ferences between the two photometric systems, which could 
significantly reduce the final T2oo/Av uncertainties (see 
Sect. 4.5.2). 



4-5.2 Effect of transformation between the ISOPHOT and 
COBE/DIRBE surface brightness photometric 
systems 

We tested the effect of the correction of ISOPHOT surface 
brightness values by the ISOPHOT-DIRBE transformation 
coefficients found in Appendix A (see Table lA*T) . Both the 
dust temperatures and the 1200/Av ratios have been recal- 
culated by the corrected surface brightness values in order 
to obtain the corrected T2oo/Av emissivities. The results are 
shown in Fig. 5b. The effect of this correction is small for 
the T2oo/Av values obtained from measurements with the 
90 and 120 ^m ISOPHOT filters. The 100 ^in-related data 
points are more displaced, due to the relatively large dif- 
ference in the scaling factors. However, the overall shape of 
the distribution did not change significantly, and the effect 
is much smaller than that of the general calibration errors 
presented in Sect. 4.5.1. 



1 




10 



4.6 The impact of various error sources on the 
emissivity — dust temperature relationship 

The distribution of data points in Fig. |5p indicates a clear 
trend in the T2oo/Av versus T relationship: T2oo/Av is en- 
hanced for T < 14 K and decreases towards higher temper- 
atures in the range 14<T<16K. In the preceeding sub- 
sections we investigated possible mechanisms which could 
explain this observation, our findings are summarized be- 
low: 

• 2MASS data provide reliable extinction maps for our 
regions. Contamination by foreground stars - which is most 
critical towards the densest regions - cannot be excluded, 
but probably has a minor effect. 

• The dense clouds in our sample, presumably also the 
cold ones, may have a ratio of total over selective extinction 
Kv higher than the standard value of 3.1. As a consequence 
Ay may decrease at most by ~19%, and the correspond- 
ing I200/AV and T2oo/Ay values increase accordingly by the 
same fraction. However, high Ry should be restriced to small 
areas in our target fields. These regions are probably ex- 
cluded anyhow due to their too low star count or too high 
extinction values (Ay > 8 mag). Therefore the likely effect of 
Kv on the final emissivity values is below ~10%. 

• Assuming a temperature dependent emissivity law (5{T) 
instead of a fixed /3 = 2 has a minor effect on the T200 / Ay 
distribution (Fig. |HJ . The differences with /3 — 2 are more 
prominent at higher temperatures, but the main trend for 
clouds colder than the dust in the DISM is not changed. 

• The possible presence of warm dust causes an overesti- 
mate of the dust temperature and, consequently, an under- 
estimate of T2oo/Ay by 10-30%. This variation is within the 
typical measurement uncertainties of the data points and 
does not alter the trend as observed in Fig. |5] 

• Calibration errors can significantly modify the T200 / Ay 
and dust temperature values. The application of ISOPHOT 



Figure 6. The histogram with shaded area shows the rela- 
tive frequency (dN/de) of the observed T2oo/Ay values rela- 
tive to that of the DISM (e = T2oo/Av/[t2oo/Av]dism, where 
[ t 20o/Av]dism = 2.4-10 -4 mag -1 ), found for our targets. Hor- 
izontal bars represent the predicted range of enhancement in 
emissivity for different types of grains following Stognienko et 
al. (1995), from bottom to top: (1) ice mantles, (2) silicate 
particle-cluster aggregates, (3) silicate cluster-cluster aggregates, 
(4) mixed (silicate and carbon) particle-cluster aggregates, (5) 
mixed cluster-cluster aggregates, (6) carbon particle-cluster ag- 
gregates (7) carbon cluster-cluster aggregates. 



surface brightness values corrected for the systematic differ- 
ences of the ISOPHOT and COBE/DIRBE surface bright- 
ness photometric systems does not change the observed dis- 
tribution of T2oo/Ay values significantly. Applying a ±9% 
general calibration uncertainty would cause a much severe 
effect, but the corrected results remain in the ~25% range 
to the original values. Transformation of the ISOPHOT sur- 
face brightness values to the DIRBE photometric system has 
only a minor effect on the Ta vs. T2oo/Ay relationship. 

From the uncertainties due to different error sources we 
estimated a general ~30% uncertainty for our T2oo/Ay val- 
ues. To match the observations by assuming a constant emis- 
sivity (that of the DISM) a 50. ..80% uncertainty is required, 
which is unrealistic, according to our investigations. The 
observed temperature versus T2oo/Ay relationship, as pre- 
sented in Fig.|2p, most probably indicates physical changes 
and is not due to artifacts in the data or to incorrect initial 
assumptions in the calculations. 
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5 DISCUSSION 

5.1 Observed variations in T2oo/Av 

We found notable differences between the original del Burgo 
et al. (2003) T2oo/Av values and the emissivities of the same 
fields reprocessed in this paper. The main differences be- 
tween our and the del Burgo et al. (2003) processing are 
that in this study we applied extinction values derived from 
2MASS data and that we used different short wavelength 
data for the temperature calculation, namely 90, 100 or 
120 nm instead of 150 /im (these new results are also in- 
cluded in Table 1 and in Fig. |2j. The dust temperatures we 
found are consistent with those in del Burgo et al. (2003). 
However, the resulting I200/AV ratios are different, and this 
discrepancy lead to notably lower T2oo/Av values (see Ta- 
ble 1). This also indicates, that it is mainly the application 
of 2MASS data that is responsible for the differences seen 
in the emissivities. 

Our values obtained for the TMC 2 (G173.9-15.7 and 
G174.3-15.9) and for LDN 1780 (G359. 1+36.7) are very sim- 
ilar to the average values presented for these clouds by del 
Burgo & Laureijs (2005) and Ridderstad et al. (2006), re- 
spectively. In these papers I200/AV and T2oo/Av values were 
derived by taking the brightness ratios from the maps in- 
stead of correlation plots. Despite the general argeement 
with our data, the detailed analysis of LDN 1780 show that 
T2oo/Av can range from 10 4 mag _1 to 4-10 4 mag _1 within 
the same cloud (Ridderstad et al., 2006). An even higher 
increase of T2oo/Av is observed in the Taurus 2 molecular 
cloud (del Burgo & Laureijs, 2005). 

An important outcome of our analysis of a large sample 
of cloud regions is that we find a clear trend of T200 / Ay ver- 
sus Td, but that the variation in T2oo/Ay is less pronounced 
than initially reported by other studies. Only three fields 
exhibit T2oo/Ay which is more than threefold the value for 
the DISM (t 2 oo/Avdism). All other regions with T d <14 K 
(12 regions) have T2oo/Ay > T2oo/Avdism- In the majority 
of the regions with Td >14K the emissivities are below 
t"20o/Avdism (10 out of 12 regions). One extreme case is 
G142.0+38.5, where 3r 2 oo/A v « t 2 oo/Avdism. 

The values of T2oo/Av in our sample refer to the aver- 
age properties (Td, and 1200/Av) over a relatively extended 
region of -~100 arcmin 2 . The enhancement observed by us of 
e < 2 in most cases suggests, that higher e values usually 
occur only in smaller regions inside the clouds, possibly in 
the densest cores. 

We have not seen large deviations in the scatter plots, 
neither in the I200/AV, nor in the dust temperature deriva- 
tion, as reflected in the uncertanties in Td and I200/AV. 
Extremely high, local emissivities are expected to show up 
in these scatter plots, and are therefore excludeable in our 
fields at the resolution we adopted, namely 3'.5. 



5.2 Changes of dust properties at low 
temperature 

Below 14 K all T2oo/Ay values are higher than 
2.4-10 -4 mag -1 , the representative value of the DISM. 
This finding agrees with the results of several earlier papers 
(see Sect. 1 and the caption of Fig. 2). The majority of our 
data points are within a well-defined range, l-to-2 times 
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Figure 7. In this figure we present model fits of T d versus 
T 20o/Ay with a dust emission of bimodal temperature distribu- 
tion. The gray curves represent different emissivities and cold 
temperatures (T c ). To derive the "virtual" dust temperature, 
surface brightness at 100 fim (short) and 200 (im (long) wave- 
lengths were used. The warm dust temperature has been set 
to T W = 17.5K. Other short wavelengths (90 or 120 fi) produce 
curves very close to the 100 ^m ones. I. e = 2.0, T C = 13.0K; 
II. e = 2.0, T C = 12.0K; III. e = 2.0, T C = 11.0K; IV. e = 1.0, 
T C = 13.0K; V. e = 1.0, T C = 12.0K; VI. e = 1.0, T C = 11.0K; 
The black dots mark our measured T2oo/Ay values. 



the emissivity of the DISM. Enhanced emissivity values 
are usually interpreted as the change of dust properties: 
coagulation of dust particles or formation of ice mantles 
on grain surfaces (see e.g. Dwek 1997). Dust grains may 
show a large variety both in composition (silicate, carbon, 
or mixed) and in structure (ice mantle, particle-cluster 
aggregates, cluster-cluster aggregates) and the different 
types show different enhancement of the emissivity. Since 
the emissivity enhancement typically is not higher than a 
factor of 2, our results are consistent with ice mantles or 
cluster of silicate particles (CCA or PCA, see Fig. On 
the other hand far-infrared emissivities produced by grains 
containting carbon aggregates are too high for most of our 
observed emissivity values. However, in some particular 
regions the observed T2oo/Av exceeds significantly the 
representative values of the DISM, and may indicate the 
existence of carbonaceous grains. 

5.3 Emissivities lower than in the DISM 

As presented in Sect. 3.3 and Fig. we find low emissiv- 
ity values for 14 K < Td < 16 K colour temperatures. Similar 
results were obtained for some specific regions by Lehtinen 
et al. (2004) and by Rawlings et al. (2005). Even though 
one can explain an emissivity lower than that of the DISM 
in terms of changed grain properties, it is also possible to 
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explain this behaviour using a bimodal temperature model 
(Cambresy et al., 2001; del Burgo et al., 2003). In this model 
we assume, that the observed far-infrared surface brightness 
is due to two emission components with different dust tem- 
peratures, both due to big grains. The observed infrared 
emission can be written as: 

If = t v>c B v (T c ) + t U)VJ B v (T v ) (3) 

where T c and T w are the temperatures of the cold and warm 
dust components, respectively, and t„ jC and t v ,w are the in- 
frared optical depths of the cold and warm components, re- 
spectively. The emissivity of the cold component may differ 
from that of the warm component, and this change is char- 
acterized by the e parameter: 

t v>c = eXT„ t tot (4) 

Tv,w = (1 — X)Tu,tot (5) 

where T Vt tot is the total effective optical depth if both compo- 
nents had identical emission properties and X is the fraction 
of the optical depth of the cold component with respect to 
the effective total opacity (for more details see del Burgo et 
al. 2003). In del Burgo et al. (2003) the cold temperature 
T c was choosen to be 13.5 K, representing the lowest colour 
temperatures obtained from their sample. However, in our 
sample we included denser molecular regions and the pres- 
ence of colour temperatures of Td « 12 K obviously requires 
a lower T c value. 

We derived the colour temperature from the slope of 
the scatter plot of two IR-wavelength surface brightness dis- 
tributions. It is easy to see that in a bimodal temperature 
distribution the variation of X causes a change in the ratio of 
surface brightness values measured at two wavelengths (see 
Fig. 8 in del Burgo et al., 2003), i.e. an effect very similar to 
the change of the colour temperature. A specific X value can 
be directly linked to a "virtual" colour temperature, as long 
as T c and T w are fixed. As we demonstrate below, mixing 
rather than change of dust grain properties can be the rea- 
son that data points in Fig. fall below the representative 
value of the DISM. In Fig.Qwe plot T2oo/Av model curves 
with mixing ratios changing in the < X < 1 range, versus 
the corresponding virtual colour temperature. In this model 
we used r„ t tot = 2.4-10 -4 mag -1 , which is the representative 
value of the DISM (Schegel et al., 1998). In all cases we kept 
T w = 17.5 K and assumed different values of T c . For e either 
e = 1 or e = 2 was adopted. 

It is clear from Fig. Q that in order to reproduce the 
higher T200 / Av values for cold regions (T < 14K) the emis- 
sivity of the cold component in the bimodal model has to 
be enhanced. 

It is noteworthy, however, that the observed low 
T2oo/Av values (14K < T < 16K), can be well fitted by using 
e — 1 for both the warm and cold components. 

This result suggests that the presence of colder com- 
ponent with the same grain properties as for the DISM is 
sufficient to explain the low T2oo/Av in a region. The most 
plausible reason is that radiative transfer in these regions 
can cause some parts to become colder while our temper- 
ature measurement is biased towards the warmer regions. 
However, a necessary condition is that the colder compo- 
nent is well mixed and not resolved by ISOPHOT, other- 
wise the scatter diagrams would become highly non-linear. 



Whether this can be achieved in a cloud model, needs to be 
investigated. 

5.4 Smaller grains 

An alternative explanation for low emissivity clouds with 
low column density is an enhancement of the relative con- 
tribution of small grains, since for a grain size a, Q a b s on a 
in the far-infrared. If this is a likely mechanism then the re- 
gions with T2oo/Av lower than the DISM would have under- 
gone some recent processing such that the relative amount 
of small grains is larger than in the standard grain size dis- 
tribution. We note that a fraction of the clouds were selected 
from the IRAS data on the basis of the high 100 /im bright- 
ness. These clouds also exhibit a high brightness in the IRAS 
12, 25, and 60 /J,m bands and could indicate a bias towards 
clouds with a grain size distibution favouring the smaller 
sizes. 



6 CONCLUSIONS 

In this paper we have analysed the FIR emission proper- 
ties in a large sample of interstellar clouds, observed with 
ISOPHOT. We have derived far-infrared emissivity relative 
to the visual extinction, T2oo/Av, for each region. The de- 
rived values of T200 / Av represent the average emissivity over 
a region typically in the order of lOOarcmin 2 . 

The derived FIR emissivities T2oo/Av show an en- 
hancement for the coldest and densest regions where 
12 K < Td < 14 K, which is most probably due to the growth 
of dust grains. We confirm a similar trend found in earlier 
papers. The enhancement of T2oo/Av for 12K<Td < 14 K 
is for the majority of the regions less than 2. This is lower 
than previously reported for some specific clouds. Our find- 
ings support models where the enhancement in emissivity is 
attributed to ice mantle growth and the presence of silicate 
agregrates on the spatial scales we investigated. 

For 14 K < Td < 17.5 K, we observe a majority of regions 
where T2oo/Av is lower than that of the diffuse ISM. FIR 
emissivities lower than that of the DISM may be explained 
by the common effect of a constant emissivity and the pres- 
ence of multiple dust temperatures along the line of sight. 
To sufficiently fit the observed values the cold temperatures 
(T c ) had to be set to ~12K without significantly altering 
the dust properties. Alternatively, these clouds in our sam- 
ple could be biased towards regions with a more significant 
small grain component in their dust size distribution. 
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APPENDIX A: COMPARISON OF THE 
ISOPHOT AND COBE/DIRBE SURFACE 
BRIGHTNESS PHOTOMETRIC SYSTEMS 

The calibrational accuracy of the ISOPHOT surface bright- 
ness photometric system is crucial for the interpretation of 
the final results. Due to the lack of proper extended standard 
objects on the sky, especially at far-infrared wavelengths, the 
only practical possibility to test the accuracy is a comparison 
of ISOPHOT data with that of COBE/DIRBE values. Here 
first we check the general relationship of the ISOPHOT and 
COBE/DIRBE surface brightness systems via the compari- 
son of P22 mini-map highly processed data product (HPDP) 
background values from the ISO Data Archive. In a sec- 
ond approach average surface brightness values of our target 
fields of the present study, obtained from ISOHPOT obser- 
vations is compared with those derived from COBE/DIRBE 
surface brightness data. 

Al Description of the COBE/DIRBE background 
determination 

The monochromatic COBE/DIRBE surface brightness value 
at a specific wavelength, date and sky position is derived 
with the PredictDIRBE tool, written in IDL 2 and devel- 
oped in Konkoly Observatory, using some routines written at 
Max-Planck-Institut fur Astronomie, Heidelberg. The Pre- 
dictDIRBE routine uses the CGIS 3 software library, devel- 
oped to help the analysis of COBE data. PredictDIRBE 
uses two main DIRBE data products, the DIRBE sky and 
zodi atlases (DSZA) and the zodi-subtracted mission aver- 
age maps (ZSMA). The detailed description of these DIRBE 
data products can be found on the COBE/DIRBE home- 
page of IPAC 4 . The main steps of the COBE/DIRBE sur- 
face brightness determination in the PredictDIRBE routine 
are the following: 

• Extraction of DIRBE ZSMA (zodiacal light 
component removed) surface brightness values for 
the 10 DIRBE photometric bands. The number of 
DIRBE pixels considered here as well as for the determina- 
tion of the zodiacal contribution depends on the extension of 
ISOPHOT maps. However, the spatial sampling of DIRBE 
measurements at long wavelenghts is very fine relative to the 
physical resolution (~30'or worse). The operations below are 
done for the median values of the DIRBE pixels taken. 

• Colour correction of the cirrus component. In 
this step the temperature of the cirrus component is der- 
ermined by fitting the 100, 140 and 240 ZSMA surface 
brightness values with a modified black-body spectral en- 
ergy distribution (SED), L. oc u B v (T), where B„(T) is the 
Planck-function at temperature T and frequency v, and (5 
is the spectral index. The 140 /im DIRBE band has a lower 
weight in the fitting process, due to its well-known noisy be- 
haviour compared to the other bands. For our calculations 
the spectral index has been fixed to (3 = 2. For A < 100 fim 
the SED is approximated by a function fitted by spline 
interpolation to the measured log(A) -log(I° IRBE ) values. 

2 Research Systems Inc., Versions 5.x and 6.x 

3 http://lambda.gsfc.nasa.gov/product/cobe/cgis.cfm 

4 http://lambda.gsfc.nasa.gov/product/cobe/dirbe_overview.cfm 



The colour correction is performed using the fitted SEDs, 
and the transmission curves of the DIRBE filters. The fi- 
nal, monochromatic surface brightness values are reached 
by the repetition of this process, until the convergence cri- 
terium is matched. We have to note, that the ZSMA sur- 
face brightness contains at least two main components: the 
cirrus (or, in general Galactic interstellar matter) emission 
and the extragalactic background. Although the two com- 
ponents have different SEDs, the difference between colour 
corrections needed for the two different SEDs is small, and in 
most cases the absolute level of the extragalactic background 
component is much below that of the cirrus. Therefore here 
we do not consider these two components separately. 

• Determination of the non-zodiacal DIRBE sur- 
face brightness at the ISOPHOT wavelength. The 
monochromatic DIRBE surface brightness values are inter- 
polated to the nominal ISOPHOT filter wavelength. For 
A < 100 /im this is performed by the spline interpolation of 
the log(A) -log(I° IRBE ) values, while for A > 100 fim the pre- 
viously mentioned i/^B^T) is applied. 

• The zodiacal component is determined from the 
DIRBE zodi atlases (DSZA products). From the date of 
the observation the actual solar elongation angle of the re- 
quested sky coordinate is calculated. In case no DIRBE ob- 
servation was performed at that solar elongation and eclip- 
tic latitude (due to the limited lifetime of the instrument) , 
the actual solar elongation is "mirrored", so that the so- 
lar elongation with the same absolute value is taken on the 
other side of the Sun, and the zodiacal component is ex- 
tracted at the corresponding dates and coordinates. This 
latter position was in almost all cases observed by DIRBE 
(for the details of COBE/DIRBE observations, see Hauser 
et al., 1998b). The zodiacal component is extracted for all 
the 10 DIRBE photometric bands. 

• Colour correction of the zodiacal component is 
performed assuming a pure black-body SED around the 
peak of the zodiacal emission, and using spline interpolation 
for notably shorter and longer wavelengths, applying other- 
wise the same iterative process as for the colour correction 
of the non-zodi component above. The resulting monochro- 
matic zodiacal surface brightness values are used to interpo- 
late to the requested ISOPHOT wavelength. 

• The final PredictDIRBE surface brightness is the 
sum of the monochromatic zodiacal and the non-zodiacal 
(cirrus + extragalactic background) components. 

A2 Minimaps HPDPs 

Mini-map mode was a 'submode' of the P22 astronomi- 
cal observing template (AOT, see Laureijs et al. 2003) of 
ISOPHOT. The mode was especially designed to obtain ac- 
curate fluxes of point/compact sources. However, apart from 
the point source flux determination, the mini-maps provide 
accurate surface brightness values of the average background 
around the target, due to the redundant observation of the 
same sky position by many ISOPHOT pixels. The mini-map 
backgrounds are assumed to be homogeneous for the whole 
mini-map area. The typical spatial extension of a minimap 
is ~4' for both the C100 and the C200 detector arrays. 

The background values can be compared with the values 
obtained by the PredictDIRBE routine for the specific cen- 
tral coordinates of the source (assuming a 30' aperture) and 
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Figure Al. Comparison of minimap background surface brightness values and background estimates based on COBE/DIRBE data at 
different wavelengths (a) 90 fim, (b) 100 jim, (c) 120 (im and (d) 200 (im. The solid lines are fitted linears, corresponding to the values in 
Table EH 



for the ISOPHOT wavelengths. Moor et al. (2003, 2004a, 
2004b and 2005) performed a full re-evaluation of a large 
sample of ISOPHOT mini-maps, containing the observations 
of normal stars, evolved stars, extragalactic and miscella- 
neous objects. The publicly available ISO Data Archive 5 
contains the re-evaluated results of these observation as 
Highly Processed Data Products (Salama et al., 2004). 

The results of the COBE/DIRBE and mini-map back- 
ground comparisons are presented in Fig. IA1I for the four 
ISOPHOT filters used in our study (90, 100, 120 and 
200 fim) . In Table IA1I we present the coefficients found for 
the correlation between the DIRBE and ISOPHOT surface 



http: / /www. iso.vilspa.esa.es/ida/index. html 



brightness photometric systems, assuming the linear rela- 
tionship 



1\ 



Sx/? IRBE + Offset 



(Al) 



where i^ SOPHOT anc j jDirbe arg ^ e surface brightness values 
at the wavelength A, measured with ISOPHOT and DIRBE, 
respectively. 

The last two columns of Table IA1I contain the mean 
relative deviations of orignal ISOPHOT and DIRBE surface 
brightness values (column #4) and the mean relative devia- 
tions of the original ISOPHOT and DIRBE surface bright- 
ness values after the correction with the fitted line ('residual 
scatter', column #5). These values are representative for the 
general relative accuracy of the two photometric systems. 

The DIRBE surface brightness calibration (as well as 
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A S Offset Ri:i Rf 

(/im) (MJysr" 1 ) (%) (%) 



90 


0.98(±0.02) 


1.87(±0.13) 


14.2 


8.5 


100 


1.09(±0.03) 


0.43(±0.13) 


16.7 


10.5 


120 


0.97(±0.06) 


0.02(±0.28) 


12.8 


11.4 


200 


0.95(±0.03) 


-0.44(±0.24) 


14.9 


10.8 



Table Al. Parameters describing the relation between the 
DIRBE and ISOPHOT surface brightness photometric systems, 
derived from the mini-map HPDP sample, assuming a linear rela- 
tionship. The columns of the table are: (1) ISOPHOT filter nomi- 
nal wavelength; (2) fitted scaling factor (S), assuming Eq. lAll (3) 
fitted offset, assuming Ea.lAll (4) mean relative deviation of orig- 
inal ISOPHOT and DIRBE surface brightness values; (5) mean 
relative deviation of ISOPHOT and DIRBE surface brightness 
values after the correction with the fitted line; 

that of ISOPHOT) contains intrinsic uncertainties, both in 
detector gain and offset. These uncertainties contribute sig- 
nificantly to the observed relative deviations listed in Ta- 
blc lAll For those DIRBE filters that are important for the 
comparison with the ISOPHOT surface brightness calibra- 
tion in this paper, the typical gain uncertainties are 13.5%, 
10.6% and 11.6% for the 100, 140 and 240 /im DIRBE bands, 
respectively (see Arendt et al., 1998). These DIRBE gain 
uncertainties indicate the presence of an additional uncer- 
tainty in the order of ~10% or below, which should account 
for the uncertainties in the PredictDIRBE estimate calcula- 
tion as well as in the intrinsic uncertainties of the ISOPHOT 
calibration. 

The typical offset uncertainties are 0.3, 0.6 and 
0.4 MJysr" 1 for the 100, 140 and 240 /jm DIRBE bands per 
pixel, respectively (Arendt et al., 1998). These are in the 
order of the offset values found for the ISOPHOT-DIRBE 
comparison (see Table IaTIl except for the 90 /im ISOPHOT 
filter, where this offset is significantly higher. However, the 
90 /jm surface brightness values may also contain an emis- 
sion contribution from a grain population different from the 
one responsible for the A > 100 /im part of the background 
SED. This effect is not accounted for in the PredictDIRBE 
routine. 

In the case of ISOPHOT, non-unity scaling factors may 
arise from imperfect knowledge of the effective solid angles 
at different wavelengths and/or from a flux-dependent non- 
linearity of the detector. ISOPHOT solid angles have been 
re-assessed recently (Abraham et al. 2005, ISOPHOT inter- 
nal calibration report) and the new values agree within 5% 
with the previous ones. Non-linearity can occur due to cali- 
bration uncertainties in the extrapolation to weak flux lev- 
els because the ISOPHOT calibration targets were always 
much brighter than the ISM in the ISOPHOT beams. The 
importance of the non-linearity can be estimated via com- 
parison with the COBE/DIRBE calibration. According to 
the results presented in Table IaT1 the relationships between 
the ISOPHOT and COBE/DIRBE surface brightness values 
are sufficiently linear. 

A3 Fields of the present study 

The target fields of the present study were observed by the 
same AOT (P22) as the mini-maps, with some notable dif- 
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Figure A2. Comparison of the average surface brightness of our 
target fields and the DIRBE surface brightness derived for the 
same sky positions and observational dates, with the Predict- 
DIRBE routine, (a) Short wavelength niters (ISOPHOT 90, 100 
and 120 /im, marked by filled circles, triangles and squares, re- 
spectively); (b) ISOPHOT 200 /im filter. The solid lines mark 
the 1:1 line. 

ferences between the two samples. First, the extension of 
the target field maps is in the 5'... 40' range, while the mini- 
maps are smaller than these. Second, the intensity of the 
target maps contains the contribution of point sources and 
also that of larger scale structures, while the point source 
emission was subtracted in the derivation of the mini-map 
backgrounds. 

The average surface brightness of our target fields were 
derived as simple average of all pixel values, for both the 
short and the long wavelength maps of a specific field. Pre- 
dictDIRBE estimates of the DIRBE surface brightness were 
calculated for the central positions of the ISOPHOT maps. 
The observational date of the maps were also considered to 
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Table A2. The 


columns of 


the 


table are: 


(1) the 


name o 



the field (central galactic coordinates); (2) central wavelengths 
of the short/long ISOPHOT filters [fim//im]; (3) short wave- 
length DIRBE surface brightness [MJysr - (4) short wave- 
length ISOPHOT surface brightness [MJysr -1 ]; (5) long wave- 
length DIRBE surface brightness [MJysr -1 ]; (6) long wavelength 
ISOPHOT surface brightness [MJysr -1 ] 



brightness scatter plots and in the determination of the 
I200/AV ratio, which is completely insensitive for offsets. 

• Not taking into account the systematic differences, 
the mean relative deviations between the ISOPHOT the 
COBE/DIRBE surface brightness values are within — 15% 
for the four investigated wavelengths. 

• The comparison of the ISOPHOT vs. DIRBE relative 
deviations with the DIRBE intrinsic detector gain uncertain- 
ties and the ISOPHOT vs. DIRBE offsets with the DIRBE 
intrinsic offset uncertainties show, that the uncertainties in 
the DIRBE surface brightness calibration have a siginificant 
impact on the ISOPHOT-DIRBE comparison. The DIRBE 
detector calibration uncertainties are at a similar level than 
other error sources (e.g. ISOPHOT detector calibration un- 
certainties, DIRBE background estimate uncertainties by 
the PredictDIRBE routine) and they may be dominant in 
some cases. 

• A representative ISOPHOT calibration uncertainty can 
be estimated by subtracting the DIRBE detector gain un- 
certainties (in average —12% for the 100, 140 and 240 /im 
DIRBE bands) from the mean relative deviations between 
the two photometric systems (in average — 15%, see column 
#4 in Table Al). The resulting value which we adopt for 
both short and long ISOPHOT wavelengths is 9%. 

• The DIRBE and ISOPHOT surface brightness values in 
the target fields show a behaviour very similar to that of the 
mini-map backgrounds. The general agreement between the 
ISOPHOT and DIRBE surface brightness values are very 
good. The relatively large deviations in some cases may be 
explained by the beamsize of DIRBE, which is large, even 
compared to the size of the ISOPHOT target maps. 



account for the annual variation of the zodiacal light compo- 
nent. The results are presented in Fig. IA2l and in Table fA~2l 



A4 Summary 

The results of the DIRBE and ISOPHOT surface brightness 
comparisons can be summarized as follows: 

• In general a very good linear relationship was found 
between the ISOPHOT and DIRBE surface brightness pho- 
tometric systems. 

• The scaling factors between the DIRBE and the 
ISOPHOT surface brightness photometric systems are close 
to unity, within the uncertainties. Even largest difference 
(100 /im) is within 9% to the DIRBE system. These scaling 
factors indeed have an impact on the derived dust tempera- 
tures and I200/AV ratios. This effect is widely discussed for 
the typical deivations in the main text (see Sect. 4.5). 

• The offsets found between the DIRBE and ISOPHOT 
surface brightness calibration based on the mini-map 
database are usually small (see Table IaTI . The largest offset 
value (-1.9 MJysr" 1 ) is found for the ISOPHOT 90 /im fil- 
ter. However, in this paper we applied the method of slope 
fitting in the derivation of dust temperature from surface 



